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Abstract 

A theoretical approach to the investigation of spin-dependent structure func- 
tions in deep inelastic scattering of polarized leptons off polarized nuclei, 
based on the effective meson-nucleon theory and operator product expansion 
method, is proposed and applied to deuteron and 5 He. The explicit forms of 
the moments of the deuteron and 3 He spin-dependent structure functions are 



*On leave from Bogoliubov's Laboratory of Physics, JINR, Dubna 141980, Russia. 
^On leave from Theoretical Physics Institute, University of Alberta, Edmonton, Alberta T6G 
2J1, and TRIUMF, 4004 Wesbrook Mall, Vancouver, B.C. V6T 2A3, Canada 
•f-On leave from Far-Eastern State University, Vladivostok, Russia. 



1 



found and numerical estimates of the influence of nuclear structure effects are 
presented. 
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I. INTRODUCTION 



The large body of experiments on deep inelastic lepton-nucleon scattering performed in 
the 70-th and 80-th, which resulted in the amazing agreement of the data with the quark- 
parton model and the confirmation of the phenomenon of logarithmic violation of Bjorken 
scaling, provided the basis for the foundation of Quantum Chromodynamics (QCD) as the 
theory of strong interactions. However, the latest experimental results in this field seem 
to display deviations from the predictions of the quark-parton model ("free QCD") and 
perturbative QCD. ( M, H). In this context, it should be stressed that basic theorems 
(sum rules) of QCD require the knowledge of the spin structure function (SSF) g™(x) of the 
neutron. 

This has been recently extracted by deep inelastic scattering of polarized leptons off 
polarized nuclear targets, e.g. deuteron (Spin Muon Collaboration (SMC) ||) and 3 He 
(E142 experiment |4]]) and these data, combined with earlier data of the European Muon 
Collaboration (EMC) on the proton [|l], are being currently used to check the predictions of 
QCD 0,0. 

The Bjorken sum rule has been computed from E142 data and, using the result 

^g{{x)dx = 0.126 ± 0.010 ± 0.015 obtained by the EMC collaboration |IJ, the value 

i 

/ (gf — g™) dx = 0.148 ±0.022. H has been found, to be compared with the theoretical pre- 
o 

diction of 0.187±0.004. At the same time, the SMC Collaboration estimated for the first time 

i 

the first moment of SSF of the "isoscalar" nucleon : M n = J dx = 0.023 ± 0.02 ± 0.015, 

o 

and combining again these data with the EMC proton data, for the Bjorken sum rule was 

found to be J (gf - <) dx = 0.20 ± 0.05 ± 0.04. 

o 

New experiments are planned at DESY J7|, CERN || and SLAC ||], also aimed at an 
improved measurement of this fundamental prediction of QCD. 

It is worth stressing here that, all information on the neutron SSF have been and will be 
obtained by analyzing DIS off polarized nuclear targets, in particular 2 H and 3 He. Therefore 
such an information can in principle depend upon nuclear structure effects, which however 
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are considered to provide only a minor correction. 

As a matter of facts, the possibility to measure the near free isoscalar nucleon structure 
functions by using polarized 2 if target is motivated by the fact that typical nuclear effects 
in the deuteron are small and are predominantly determined by well-known spin-orbital 
structure of the deuteron wave function (c.f. Ref. H§]). At the same time, the interest in 



polarized 3 He targets stems from the observation [10] that since the dominant component 
of the 3 He wave function is fully symmetric S wave with the two protons in a spin singlet 
state, a polarized 3 He can be associated to a large extent with a polarized neutron. It should 
however be pointed out that the precision required by a significant check of the Bjorken Sum 
Rule would demand a quantitative estimate of all possible nuclear effects. For example, in 
case of 3 He, the S' and .D-wave components of the small realistic three body wave function 
generate a proton contribution to the 3 He polarization and asymmetry, which has to be 
subtracted in order to obtain information on the neutron properties. Recently, the proton 
contribution to the process 3 He(e, e')X has been quantitatively evaluated by full convolution 
approach, where not only Fermi motion [JTlj] but binding as well were taken into account 
by generalizing the usual convolution approach (see e.g. Ref. [0), by introducing the 
concept of spin dependent spectral function jrj,[K|], and nuclear effects have been indeed 
found to be small. Relativistic light cone calculations have been performed in Ref. fl6 



using the spin dependent spectral function of Ref. 15]; the obtained results practically do 



not differ from the non relativistic ones obtained in Ref. |L4 . 

The underlying reason for the applicability of the convolution model is the existence of 
two typical momentum scales in deep inelastic processes, which leads to the factorization of 
the amplitude of the reaction into two pieces, depending respectively on the "large" external 
momentum, i.e. the structure function of the nucleon, and the "low" typical momenta of 



nucleons in the nucleus, i.e. the momentum distribution of the nucleons ||17|| . A rigorous 
method to analyze such a factorization is based upon the Wilson operator product expansion 
(OPE). A theoretical approach to investigate deep inelastic scattering on nuclei by using the 
OPE method within the effective meson-nucleon theory with one-boson-exchange (OBE) 



interaction has been suggested in Ref . [|18[|19f ; these calculations have been performed within 
well-defined approximations, such as the leading twist approximation in the OPE and the 
lowest order approximation on the meson-nucleon coupling constant. Such an approach 
allows one to derive a convolution formula which includes binding effects and, at the same 
time, preserves the energy-conservation sum rule ||12|| . 

In the present paper the model is extended to deep inelastic scattering of polarized leptons 
off polarized nuclear targets. To this end, the set of operators, providing the basis for OPE, 
is extended by considering the axial operators in terms of nucleon fields interacting with 
meson fields. Using the non-relativistic reduction, an explicit form of the operators relevant 
to describe polarized deep inelastic scattering are found. A particular attention is paid 
to the investigation of the properties of the physical nucleon, e.g. the physical mass, the 
meson cloud, the renormalization constants and the SSF's which appear in the calculation. 
The explicit expressions for moments of SSF of the lightest nuclei, the deuteron and 3 He, 
are found and the inverse Mellin transform reconstructs the corresponding nuclear SSF in 
the form of convolution of the nucleon SSF and effective distributions of the nucleons in 
a nucleus. The obtained formulae are generalized to the case of heavy nuclei and found 
to be similar to those used in the conventional convolution approach |T2|] . As an example 



of applications of the method, the SSF of the polarized deuteron and 3 He are numerically 
estimated and compared with recent experimental data. 

The paper is organized as follows. In Section II the basic formalism is presented. The 
antisymmetric part of the Compton scattering amplitude is defined in terms of the axial 
twist-two operators in the OPE method and their explicit form in the non-relativistic limit 
is computed. In Section III the moments of the spin-dependent structure function gi(x) for 
the physical nucleon are evaluated in terms of the corresponding moments of bare nucleons 
and meson cloud contribution. The moments and the SSF g^{x) and gl He (x) are calculated 
in Section IV, where an extension to complex nuclei is proposed and a formal comparison 
with the conventional convolution approach is illustrated. Preliminary results for A = 2 
have already been presented in Ref. |2(| . The results of numerical calculations are presented 



in Section V and the Conclusions in Sections VI, respectively. 



II. BASIC FORMALISM 



A. Kinematics and notation 



The spin-dependent structure function can be determined experimentally by measuring 
the asymmetry in the reaction with polarized particles 

l+A — ► V + X. 

In the one-photon exchange approximation the relevant part of the cross section can be 
written in terms of the antisymmetric parts of the leptonic L^ u and hadronic W^ u tensors 

^ ^ E ' L [H W^ (11.1) 



dE'dVL Q 4 M A E 

where [piu] are the antisymmetric indices, E and E' are the initial and final energies of the 
lepton, Q 2 = — q 2 is the square of four-momentum transfer, a stands for the electromagnetic 
fine structure constant, and M A is the mass of the target. The leptonic tensor is obtained 
from Quantum Electrodynamics and the antisymmetric part of the hadronic tensor W^ v is 
expressed in terms of two independent spin structure functions G± t 2 

W [H (p A , q, S) = ie, uXu q x (g^, Q 2 ) S° + [(p A ■ q) S° - (S ■ q)p° A ] ^jj^) , (H-2) 

where p A denotes the 4-momentum of the target, q stands for the time-component of the 4- 
vector q (in the rest frame of the target we use the notation q = v) and S is the polarization 
four- vector, normalized as S ■ S — — 1 and satisfying the relation S -p A = 0. In what follows 
we will consider deep inelastic scattering off polarized targets with spin one-half (the nucleon 
and 3 He) and one (the deuteron), in which case the vector S may be computed in quantum 
field theory as the mean value of the canonical spin operator defined by Noether's theorem 
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nj.[] We choose the quantization axis in the opposite direction to the photon momentum, 
g" = (i/,0,0,-|qi). 

In order to measure the SSF G\ and G2 appearing in eq. ( |11.2| ), one has to consider the 
difference between the cross sections corresponding to parallel and antiparallel electron and 
target spins, respectively, for one has 

where 9 is the electron scattering angle. In the Bjorken limit (Q 2 — > 00, v — > 00, x = 
Q 2 /{2paC[) is fixed), which will be considered from now on, the functions G12 are predicted 
to depend only upon x, yielding the "true" dimensionless spin-dependent structure functions 
gi(x) = vGi(v, Q 2 ) and gi{x) = z/ 2 G 2 (z/, Q 2 )/Ma- Defining the asymmetry as the ratio 

d 2 a/dE'dn (tt) - d 2 a/dE'dn (T|) 
ll[X> d?a/dE'dn{]]) +d 2 a/dE'dn(U) ' 1 ' ' 

one obtains 

where -^(x) is the usual spin- independent structure function. The asymmetry ( [11. 5| ) and, 
consequently, the SSF are the main objects of experimental and theoretical investigations. 

Since any structure function, either spin-independent or spin-dependent, can be repre- 
sented as a linear combinations of helicity amplitudes, it is instructive to perform further 
analysis in the helicity basis. For a thorough analysis of helicity amplitudes in case of tar- 
gets with an arbitrary spin we refer the interested reader to [p^] . We define the amplitudes 
A\m,\'M' as 

AxM.X'M^eTWwiSyy, (II.6) 



1 For example, its components are S\ = (0, 0, 0, ((t z )) for the nucleon and 5a - 
-(i/M^ex^pS^S^pp, for the deuteron with total momentum pd and polarization 4-vector S^. 
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where A, A' and A4, M! are the spin projections of the photon and target, respectively, along 
the z-axis, is the polarization vector of a helicity A photon, e± = =f(0, 1, ±i, 0)/v2, Sq = 
(— |q|, 0, 0, v)/\/Q 2 , and the other notations are self-explaining. A simple calculation 
shows that the structure function g\ in the Bjorken limit is defined by the amplitudes ( [11. 6|) 

Q 2 

A+±,+± = T9i ± ^ 92 (II. 7) 

and it reads 

<7i = -^(A ++ , ++ -A+_, + _), (II.8) 

Actually, the amplitudes ( |II.6| ) may contain the symmetric part too, which involves the 
structure functions Fi(x, Q 2 ) and, for spin one targets, &i(x, Q 2 ) P^j . In our case this part 
is irrelevant, since it does not contribute to the asymmetry ( |11.5| ) and to the SSF. 

The SSF, as well as the hadronic tensor W^], can be directly calculated, via the optical 
theorem, from the imaginary part of the amplitude for forward Compton scattering of virtual 
photons off hadronic targets 

W [H = ilmr w . (II.9) 

In what follows the Compton amplitude Tj^i for nuclear targets and its relation with SSF 
through eqs. ( |II.6|) and ( |II.8|) will be obtained. 

B. The Compton scattering amplitude for nuclear targets 

To establish a direct connection between SSF and the amplitude Tj^j we decompose the 
latter in the same form as ( [II. 2| ), viz. 

T [H = e^g* l - ( ai (x,Q 2 )S° + ((p A ■ q)S° - (S ■ q)jf A j) , (H.10) 

where the Compton spin-dependent structure functions ai^(x) are related to the Compton 
helicity amplitudes 

h\M,\'M' an d the deep inelastic SSF by the following relations: 
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h XM ,x>M> = eTT [tlu] (S)e u x , (11.11) 

oci{x) = ~(h ++>++ - /*+_,+_); (11.12) 

9i( x ) = 77- Im on{x). (11.13) 

Z7T 



Using the dispersion relation for the function ot\(x) and eq. (|II.13| ), we get the expression 
for a\(x) in terms of the moments of the SSF g\ 

oo / 1 \ n +l rl 

ai (x)=4 £ - / dyy n gi (y). (11.14) 

n=0,2,... J ° 



The integral in the r.h.s. in eq. ( [II. 14 ) is the n + 1 moment M n+ i(Q 2 ) of the SSF. We will 



use eq. ( 11.14 ) to obtain the structure functions from the explicit expressions of the nuclear 
Compton amplitude T^ v y 

The computation of the amplitude TL„i requires the treatment of two relatively inde- 
pendent questions: i) the analysis of the properties of the time ordered product of electro- 
magnetic current operators at high momentum transfer (Q 2 — > oo), which characterizes the 
short distance physics, and ii) the determination of the vectors of the nuclear ground state 
\Pa)i which essentially characterizes the large distance physics. This is seen explicitly from 



the expression of the amplitude (|II.10|) which is of the form 



7[ H = i J d^exp( iq O( P AM\T (UOUO)) \VaM). (11.15) 

The behavior of the T-product of the electromagnetic currents may be established in a 
general form directly from eq. (|II.15|) . In the Bjorken limit (Q 2 — > oo) the main contribution 
in the integral (|II.15|) comes from small space-time intervals, £ 2 — > 0. In this limit the 
arguments of electromagnetic currents coincide and T-product contains singularities. A 
consistent method to analyze these singularities is based on the Wilson operator product 



expansion P3[ . According to the OPE on the light-cone p4| , the product of two arbitrary 
operators A and B factorizes into two pieces when their arguments are separated by a small 
space-time interval; the first piece contains singularities (the c-number coefficient functions, 
or Wilson's coefficients) and the second one appears as a set of regular local operators, which 
are well-defined in field theory. Then the local operators are expanded in a series 

9 



A(0iJ(O) ~ ^C^(O). (IL16) 

n 

The operators O^ 1 "^™ are defined to be symmetric and traceless in all Lorentz indices 
fix . . . fi n . The quantity n is the Lorentz spin. Another quantity, the twist, defined as 
t — d n — n (d n - is the canonical dimension of the operator O n ) plays an important role 
in the theory of deep inelastic scattering processes. Namely, only the lowest values of r 
contribute to the matrix elements of the Compton amplitude [|17]]. Therefore in the leading 
order of the twist (r = 2), the r.h.s. of eq. ( |II.15D can be rewritten as 

where t tags the fundamental fields of the theory under consideration and oJ tT ' u '"' 4n ^(0) are 
the relevant twist-two operators constructed from this fields. The transformation properties 



oo / 2 

T [H =^aV E C n,t(Q 2 ) 75 q» 1 ---q f > n (PAM\6}' T ^ } (0)\p A M), (11.17) 



t;ra=0,2, 



of the amplitude Tj^j restrict the Lorentz spin, (n + l), in eq. (|II.17|) to take only odd 
values. 

It is worth emphasizing that in ( |IL17|) , due to the factorization in the OPE, the coefficient 
functions C n are related to short distance, "subhadronic" , physics (depending on the large 
momentum q), whereas the matrix elements of the operators cK 0- ^ 1 — characterize the large 
distance physics (depending on typical nuclear momenta). 

In order to separate the part contributing to gi(x), it is convenient to rearrange the 
symmetric Lorentz indices {crui . . . u n } so as to obtain two operators, one with no definite 
symmetry Q a {^-^n} ^ an( j foe other with mixed symmetry [^5[ 



viz . 



I n + l 



(11.18) 



Lt=l 

Then gi gets contributions from O^i^---^} on iy The mixed symmetry operators enclosed in 
the square brackets of eq. ( |11.18j ), define the twist-2 part of g2(x), the so-called Wandzura- 
Wilczek contribution [ ERCT . In this paper we consider the amplitude Tj M „] defined by 
eqs. ( |1I.10| ) and ( EI.17 ) with only the first operator from eq. (|II.18|) , that is, the part 
concerning the structure function g±(x) only. The structure function g2^x) is expected to be 
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very small, (in the parton model it is exactly zero), and will not be considered here (for the 
investigation of nuclear effects in g 2 (x) see, for instance Ref. |27|). 

The form ( |11.17| ) for the amplitude ( |11.15| ) near the light-cone is valid in the framework of 
any renormalizable field theory |23| , |2^f| . Thus, the problem of its analysis is formulated now 
as a consistent calculation of both the coefficient functions C n> t and the matrix elements of 
the twist-two operators O^ 1 '"^ sandwiched between nuclear ground state vectors \pa). So 
far there does not exist a realistic field theory by which one could compute simultaneously 
both pieces, since if one of them is calculated in a more or less self-consistent way, then the 
other one relies on phenomenological approaches. 

For example, within perturbative QCD, because of asymptotic freedom it is possible to 
analyze the properties of the Wilson's coefficients C ntt (a s , Q 2 ) by computing directly, at least 
in principle, the corresponding Feynman graphs up to the desired order in a s |28| . 



Making use of the renormalization group equations, the Q 2 — dependence of these co- 
efficients may be established as well. However, in practical calculations of the structure 
functions and comparison with experimental data one needs the matrix elements of the 
operators O t (0), which are related to the non-perturbative region of QCD and, hence, are 
parameterized from the experimental data. Since the basis of operators in OPE, eq. ( |L1.17| ), 
and the Wilson coefficients are target independent, all information about the target is con- 
tained in the unknown matrix elements. Therefore in the treatment of the "QCD motivated" 
models of nuclear effects in deep inelastic scattering, one is forced to introduce parameters, 
which neither can be computed theoretically from QCD, nor can be fixed from independent 
experiments [p9f[ . 

The problems related to the short and large distance are formally solved, in case of deep 
inelastic scattering, within non-asymptotically free theories with spinor fields (nucleons) in- 
teracting with the massive bosons (mesons) via pseudoscalar and vector couplings |30| , |3T|j . 
These field theoretical models with a renormalizable interaction allow a perturbative in- 
vestigation of the Wilson's coefficients and of the corresponding matrix elements, by sum- 
ming the leading logarithmic corrections. These examples can be considered as an idealized 
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meson-nucleon theory, since the realistic models necessarily involve phenomenological cor- 
rections, such as vertex meson-nucleon form factors, effective coupling constants and meson 
|32| , |33|| . In realistic meson-nucleon models, the exact results of Refs. pO| , |3T | 



masses 



serve as 



a hint for a formal approach and further phenomenological adjustments. It is obvious that 
the target independent Wilson's coefficients are not calculable within such a theory and they 
ought to be parameterized from experiments, for instance, from the experiments on deep 
inelastic scattering off free nucleons. 

We apply the effective meson-nucleon theory and the OPE method to deep inelastic 
scattering off nuclei (c.f. Refs. |^JT8|JT9|j27[| ) . This theory allows one to describe fairly well 
the NN-interaction at relatively small energies, the nuclear bound states \pa), the binding 
energy and other properties of light nuclei |35| . Then, the axial twist-two operator 
for the system of interacting nucleon and meson fields may be written down in the form 



O£{m-^}(0) 



x {:iV(0)7 <7 7 5 d w N(0):}, 



[11.19) 



where N stands for the nucleon spinor fields. Note the implicit presence of the meson degrees 
of freedom, via interaction, in the definition of the operator O^ A ' 1 " Mn ^(0). In order to obtain 
the explicit expressions of the operator ( |11.19| ) and calculate its matrix elements, we need 
the Hamiltonian of the system. In order to achieve self-consistency, this Hamiltonian has 
to provide simultaneously the equation of motion for the interacting fields and the target 
ground state: 



N = i[H,N], 
H\ PA ) = M A \p A ). 



(11.20) 
(11.21) 



Equation ( |11.21| ) for the nuclear ground state has been solved within the non-relativistic 
limit using effective Hamiltonians containing n, a, uj, p, rj and 5 mesons (one-boson-exchange 
(OBE) approximation) p2]fJ5|] . Consequently, the operators ( |11.19| ) are also to be calculated 
within the non-relativistic limit. The strategy of our calculation is therefore as follows: i) 
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we choose the appropriate covariant Lagrangian, giving the classical equation of motion 
for interacting meson and nucleon fields; ii) these equations are reduced non-relativistically 
and the effective non-relativistic Hamiltonian is obtained; iii) using the same procedure of 
non-relativistic reduction, we compute the explicit form of the operators ( |11.19| ). 

The procedure of non-relativistic reduction of classical equations of motion for the inter- 
acting meson and nucleon fields has been established in a number of papers and could be 
found in details, for instance, in Refs. [|3~6l ,[37ll . 

Below we perform explicitly the calculations with pseudoscalar-isovector coupling, the 
pion-nucleon interaction. The result is generalized to the case of other kinds of couplings. 
Introducing the isospin formalism we redefine the Wilson's coefficients in eq. (|II.17| ) as a 
diagonal (2 x 2) matrix in the isospin space with the proton and neutron coefficients on 
the main diagonal, (C n ) = C na 5 a 8, a, 3 = 1,2, and use them into the definition of the 
operators ( |II.19|) . 

The contribution of the operators O cr ^ 1 - ^ n ^ ( |11.19| ) to the Compton helicity amplitude 
then becomes 

h+M,+ M = - E ( Q2 J (^n)A' 



n=0,2,. 

2\ 



(0+) A = y ( PA M\ :iV(0)C n (Q 2 ) 7 + 75^ N(0): \p A M), (11.22) 

where 7+ = 70 + 7 2 and <9_ = <9 — d z . 

We perform the non-relativistic transition of the fields iV(0), following the method de- 
scribed in Refs. p6|j37[| , and, using the equation of motion ( |II.21| ), we compute their n-th 
cL derivatives. The resulting operators, being composed from interacting nucleon fields, 
explicitly involve the meson degrees of freedom. Skipping some rather cumbersome details 
of calculations we write below the explicit non-relativistic form of the operators 0+ as a 
sum of operators up to second order of growth coupling constant g n , up to g\: 

On = Nn + + N ^ n , (11.23) 

X N, n = m n I ^^JV(% 1)P2 )[£V,P 2 )] a + ( Pl , Sl )C n a(p 2 ,s 2 ), (11.24) 

J (27r) L Js i s 2 
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o 



A 

Nir,n 



m 



ig n r dpidp 2 dk 



2m J (2vr) y J 2 u(k) 



(jVl-)(k)fet(k)-iVW(k)6,(k) 



S A (pi,p 2 ) cr • k a + (pi, s^&nT 3 a(p 2 , s 2 ) + h.c. 



w NN,n 



-m 



SI S2 

2 ' dpidp2(ip <ik 



(11.25) 



4m 2 7 (2tt) 12 2a; (k) 



S A (pi,p 2 )<r -k 



>(k) + iV«(-k) 



[cr-k] : a + (pi,Si)C , n r J a(p 2 ,S2)a + (p,s 3 )r i a(p + k, s 4 ) : +/i.c. (11.26) 



where m is the mass parameter of the Lagrangian, the mass of bare nucleons, Tj are the 
isospin matrices, a ± (p, s) are creation and annihilation operators of a bare nucleon, proton 
or neutron, with spin s and momentum p, fe - (k) create or annihilate the j-th pion with 



X%. sA Xs 2 ; Xs is the 



si s 2 



momentum k; uj = \jm\ + k 2 , (m n is the pion mass) 
Pauli spinor and the operator £ is the non-relativistic analogue of the four-vector spin 
operator for the spin-^ particles 

2^cr • (pi + p 2 ) , A = 

- A ( 1 " E1 ^) + 

+ { a ' Pi P2 + a ■ P2 Pi ~ Pi • P2 ^ A - i\pi x p 2 ] A ) , \ = x,y,z 
The functions iV^ in eqs. ( |II.24|) -( [II.26 ) depend only upon momenta and are of the form 

^ 1) (Pl,P2) 



S A ( P i,p 2 ) 



(11.27) 



\ + Pi 2 + P 22 + +^2 2 



4m 2 



2m 



NW(k)=± 



[ 1± ^+( k )~ 


"-) 


2m 





w+ (k) 



1 ± 



fc 2 
2m 



(11.28) 
(11.29) 
(11.30) 



where uj + = c<j + The operators ( |II.24|) -( [[I.26| ) are the basic result of non-relativistic 
OPE method within the OBE approximation. Their matrix elements will determine the 
polarized deep inelastic scattering of leptons on those nuclear targets, e.g. the two and three 
nucleon systems, which are well described within the effective meson-nucleon theory with 
OBE potential. 

Note that our method is based upon perturbation theory within the effective meson- 
nucleon theory. Since the effective meson-nucleon coupling constants are large, it 
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can be argued whether a perturbation expansion can be applied at all. We simply follow 
here the customary use of perturbation theory in the effective meson-nucleon theory p2| , |33 



justified by the success of OBE model. The price one has to pay in such an approach is the 
use of some phenomenological ingredients. For example, since the nuclear forces are strongly 
repulsive at small distances, the physics at such distances is mocked up by a short-range 
repulsive core, which is handled partly by introducing meson-nucleon vertex formfactors. 
Thus, in the calculation of the corresponding matrix elements the uncertainties of the nuclear 
force at very short distances are insignificant. Analogously the short distance contributions 
to the nucleon and nuclear structure functions are hidden into the Wilson's coefficients, 
which are the new "effective constants" in the effective meson-nucleon theory |18| , |l9 , |27 ,34 



From the QCD point of view this situation corresponds to the picture when hadrons, 
viewed as a core of valence quarks surrounded by a sea quark-antiquark pairs and gluons, 
are approximated by a core of bare nucleons with a correlated color neutral quark-antiquark 
pairs, the meson cloud. In OPE such a picture means that much of detailed dynamics 
of the quarks is embedded in the coefficient functions C n (Q 2 ), the influence of the meson 
cloud is rather included into the matrix elements of the operators O a ^ 1 "^ n \ In case of the 
effective meson-nucleon theory a preliminary investigation of the twist behavior of eq. ( |11.17| ) 
is hindered by the fact that all the renormalization effects are included into the effective 
constants of the theory (coupling constants, meson masses, vertex formfactors etc.), so that 
the renormalization group equations are here, in a sense, inapplicable. The use of the 



twist-two as a leading term in eq. ( fLT? ) is to be regarded as an assumption and it should 
be verified "a posteriori", by looking at the quality of the final results and, possibly, by 
investigating the higher twist corrections. 

III. MOMENTS OF THE NUCLEON STRUCTURE FUNCTIONS 

As mentioned above, the non-relativistic expressions for the operators ( |11.24| )-( [11.26| ) and 
the Compton amplitude ( |11.22j ) involve, in particular, the bare (unknown) parameters, i.e. 
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the mass parameter from the Lagrangian of the theory and the Wilson's coefficients C n (Q 2 ) 
which are not calculable within our approach. Whereas the mass parameter may be fixed by 
introducing into the Lagrangian a corresponding counterterm, the coefficients C n (Q 2 ) are to 
be related with the SSF of the physical nucleons. Therefore we first investigate the ground 
state and the SSF of the physical nucleons. 



A. Nucleon ground state 

For the investigation of the nuclear ground states in the effective meson-nucleon theory it 
is convenient to use the Tamm-Dancoff method. For a physical nucleon with an isospin index 
a, momentum q and a given z-projection of the spin s, the Tamm-Dancoff decomposition is 
given by 



\N) n , a , s = V 1 " Z NVo \N)^^ s + <^i|A%) qiaiS + . . . , (III.l) 

where is the normalization constant defined by the condition (N\N) = 1 and \N) and 
\Ntt) represent the basis vectors of the states with one bare nucleon, and with one bare 
nucleon and one meson, respectively. The coefficients tpi in the expansion ( |IH.1|) are the 



operators in the momentum space and define the corresponding wave functions 

<^o|iV) q , QiS = a+(q, S )|0), (III.2) 

<pAN*:) W = / ^?^;r:' J ( P ,k,q)ai(p, S ')&;(k)|0). (IH.3) 
J (2tt) 

The explicit expressions for the wave functions <pi is found from the condition that the 



\N) is the state of the physical nucleon , which means that it obeys eq. ( 11.21 ) with the 
physical nucleon massQ , viz 



2 In eq. QHI.l ) we keep only the first two terms. The next term, if2 \ Ntttt) , in spite of being 
proportional to g 2 , does not contribute to the nucleon SSF, as it can be seen from eqs. ( 11.24] )- 

(TO )- 
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A:?(P,k,q) = -(27r)V 3 ) (p + k-q) 



[a ■ k} s , 



s' M 



(III.4) 



2m N u(k)yj2u(k) 

In particular, as it can be seen from ( |ill.4| ), in the non-relativistic approximation there is 
no interaction between the nucleon and the pion with the angular orbital momentum 1^1. 
Finally one can find the mass counterterm 5m and the renormalization constant that 
are of the form 

dk 3k 2 „ ql r dk 3k 2 



5m 



9l 



4m 2 J (2tiY 2u 2 (k)' 



4m 2 J (2%y 2uj 3 (k) ' 



(III.5) 



B. The matrix elements 



Now we are in the position to compute the matrix elements of the operators ( II.24j) and 
( |II.25| ) with nucleon ground state vectors ([III.1D . The operator (|II.26| ) is of two-body origin, 
hence it does not contribute to the nucleon matrix elements. 

Schematically the matrix elements for a nucleon at rest are as follows 



M.E. ~ a ,sAN\^Y,°N,nV0\N)a,s z (1 " Z N ) 
n 

n 



(IA + renorm. 
(recoil) 
(interaction) 



(III.6) 
(III.7) 
(III.8) 



The four different matrix elements given by eqs. (|III.6|) - (|III.8| ) are known, respectively, as 
the impulse approximation (scattering off bare constituents), the renormalization and recoil 
contributions, and the term of pure interaction origin (self energy-like correction). The 
contribution of these matrix elements with given isospin to the helicity amplitude ( pi.22| ) , 
sandwiched between the states a, may be explicitly written in the form 



2 ++ 



- (C n 



E 

n=0,2,. 

9l 



Cr 



l\ n+1 

xj 

dk 2k 2 -k 2 



5m 

m N 



«« 4m 2 N \J (27r) 3 u 2 (k) 



2a; (k) 



9l 



aa 4m 2 N 



f/k ivW(k) 3k " 



(2vr 



LU 2 (k) 
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(III.9) 
(111.10) 
(III.ll) 



where a = 1, 2 corresponds to the matrix elements on proton and neutron respectively. 

The physical meaning of the obtained result becomes clear if eqs. ( |111.9| )-( [Ill.ll| ) are 
depicted in terms of Feynman diagrams. Fig. 1 represents the helicity amplitude for the 
proton. It can be seen that the proton amplitude is determined not only by the proton Wilson 
coefficients C n> i, but by the neutron coefficients C n ^ as well multiplied by a factor 2 which 
comes from the relation between the pion-nucleon coupling constants: g n ±N = V^fiVTV- 
From eqs. QII.13| ) and (|II.14|) it is easy to show that in the present approach the Wilson's 
coefficients in the OPE are proportional to the moments of SSF of the bare nucleons (proton 
(p) on neutron (n)) 

a,i( 2 ) =2M n+1 (^ (n) ), (111.12) 

with 

i 

M n (/) = Jdxx^fix). (111.13) 



Relations ( |III.9| )- (|III.13|) determine the moments of SSF of the physical nucleons in the 



effective meson-nucleon theory. Accordingly, the moments should be parameterized in order 
to describe the scattering off free physical nucleons, and they should be replaced by this 
parameterization whenever they appear in nuclear matrix elements. The first (n = 0) 
moments of the nucleon SSF play an important role in deep inelastic scattering of polarized 
particles, since they define some integral relations among the SSF, known as sum rules. As 
n = the interaction term in eqs. ( |IH.10|) -( |IH.11|) vanishes (see also eqs. ( PX27| )-( |II.29|) ) and 



one obtains 



Mf ^) = M? (n) (l-Z 



N) 

p(n) oM n (P) 



+ (M? w + 2M? 

gl r dk 2k 2 z -k 2 
ter&J (2tt) 3 o; 2 (k) 



x ±5r / ,„ ^3 (ni.14) 



The second term in r.h.s. of eq. (Jill. 14j) is the contribution of recoil diagrams, which is 



partially canceled by the renormalization term |pq|. The remaining part is the additional 



counterterm to renormalize the composite axial operator ( |II.19| ). For n = the operator 
( |11.19| ) is the spin operator for the spinor fields and its matrix element, ( |111.14j ), is the mean 



value of the spin projection over polarized states of the physical nucleon. Equation ( |111.14| ) 



demonstrates that the spin projection of the physical nucleon differs from the mean value 
of the spin projection of the bare nucleons. This is expected, since the interaction between 
the core and the meson cloud, with orbital momentum 1 = 1, slightly redistributes the spin 
among the constituents. The total angular momentum of a meson-nucleon system, is a sum 
of the orbital momentum / = 1 and the spin s = 1/2, so that both parallel and antiparallel 
polarizations of the core contribute to the nucleon polarization. From the point of view of 
the present approach, the orbital moment of the meson cloud affects the spin distribution 
of bare nucleons inside a polarized physical nucleon. The effects of the meson cloud on 



unpolarized deep inelastic scattering on nucleon have been investigated in refs. |39].|4(| and 
are known as the Sullivan processes. 

Another interesting consequence of the eq. ( [111.14| ) may be derived by analyzing the 
difference between the first moments of SSF of the proton and neutron, the Bjorken sum 
rule. In this case one finds that the Bjorken sum rule on free nucleons may be computed by 
considering the bare constituents and the meson cloud corrections 



/ (g P i^Q 2 ) ~ 9i(x,Q 2 )) dx = M?(Q 2 ) - M"(Q 2 
= (M?(Q 2 ) - M"(Q 2 



(ni.15) 



3m 2 N J (27r) 3 cu 3 (k) J ' 

where the r.h.s. of eq. (|III.15| ) is the BSR for bare nucleons corrected by the interaction with 



the meson cloud. Explicit numerical estimates of the role of the meson cloud in polarized 
deep inelastic scattering processes on nucleons is ambiguous and requires a proper investi- 
gation. In fact one needs a consistent procedure for the regularization of integrals appearing 
in eq. ( |111.15| ) and for the inclusion of the vertex form factor into the meson-nucleon vertex. 
Also, some assumptions about the properties of the bare nucleon core are necessary. The 
role of the orbital momentum of the virtual pions in determining the spin of the nucleon has 
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to be investigated as well. These problems are beyond the goal of the present paper and 
will be considered elsewhere. 



IV. NUCLEAR MATRIX ELEMENTS 

In the previous section we explained in details the basic formalism and apply it to the 
physical nucleons. The basic idea of the calculations of nuclear matrix elements remains the 
same: once the non-relativistic expressions for the axial operators have been established in 
terms of nucleon and meson fields, the vectors of nuclear ground states are to be defined in 
the same manner, i.e., making use of the same non-relativistic Hamiltonian which has been 
used to derive eqs. ( |11.24| )-( JL11.3D . In what follows we keep the terms up to the second order 



in the meson-nucleon coupling constant g, which corresponds to the usual approximations 
in nuclear physics in deriving the potential and Schrodinger equation, g 2 - approximation. In 
this sense the present approach pretends to be self-consistent. The condition of a consistency 
is that the Hamiltonian, with OBE interaction, should indeed describe the real nucleus. We 



choose the realistic OBE potentials, such as the Bonn |35| or Reid [f|l[] ones, which give a 
good description of light nuclei. 

Below the explicit expressions for the moments of the 2 if and 3 He SSF's, are derived. 

A. The deuteron 

The vectors of the deuteron ground state in the Tamm-Dancoff approximation are given 
by a relation similar to eq. viz 

<p \NN > = ^ 1S2 (Pi, Pa) (2tt) 3 5( Pi + p 2 ) a+( Pl , Sl )^a+( P2 , s 2 ) |0 > (IV. 1) 

(pi \NNn > = 

ip? 8 *** ( Pl , p 2 , k) (2vr) 3 5( P i + p 2 + k) a+( Pl , si)a+( P2 , s 2 )6+(k) |0 > , (IV.2) 

where e a ^ is the Levi-Civita tensor which describes the isospin function of the deuteron and 
the nucleon spins s±, s 2 and orbital momentum in (Pq iS2 ( P i, P2 ) are summed to the total 
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angular momentum J = 1. The quantum numbers a, (3 and j are combined so as to give the 
total deuteron isospin T = 0. The physical meaning of the coefficient (fo can be understood 
by projecting the bare deuteron state onto the state with two nucleons located at points 
ri, r 2 in co-ordinate space 

< r 1; r 2 \^ Q N l N 2 > 

= J 7^ e * kiri+ * k2r2 ( 27r ) 35 (Pi + I») ^ (Pi> X.-xX?,^^ 
x (0|aa(ki,si)a / g(k2,S2)a^(pi,si)aJ(p2,s 2 )|0), (IV. 3) 



It is clear from eq. (|EV.3|) that y2o(p,p) with spin, Xs, an d isospin, r) a , functions is 
the conventional deuteron wave function in the momentum space, obeying the Schrodinger 
equation with OBE potential PS| . In what follows a more conventional notation for ipo will 
be adopted, viz. y?o(p,p) = ^m(p)- F° r the function ipi we have 



Pi 



Pl , P2 ,k) = -(27T)V 3 ) ( Pl + P2 + k) 



2m ^v o;(k)^/2a;(k) 

x |^ 1S (pi,Pi)— [a- ■ k] SS2 [r^ - (1 - 2) | . (IV.4) 

The renormalization constant for the deuteron state contains the nucleon contribution 
1II.5|) and the exchange part, Z D = 2Z N + Z D , where 



We proceed now with an analysis of the moments of the deuteron SSF. To begin with, 
let us define the isoscalar nucleon SSF, g^(x) = (gx(x) + g™(x))/2. Then the n-th moment 
of the deuteron solely depends on M n (gi). Using eqs. ( |111.9| )-( P"11.11| ) we get 

M„ +1 ( 9 f) = M»«tof ) {1-4^/ (i^^k) ( N » +) ( k ) k2 + ( 2k * ~ k2 ))} 

+ M„ +l(9 f ) \-Z N + (2*2 - k») + . (IV..) 

The moments of the deuteron SSF's are determined by matrix elements similar to those 
in eqs. ( [III.6|) -( [iII.8|) . In the nuclear case all the operators ( |L24D -( |L26 ) contribute to the 
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corresponding matrix elements. While evaluating the matrix elements one obtains contribu- 
tions corresponding to eqs. ( |111.6| )-( |I11.8| ) computed with the deuteron wave functions ( |LV.1|) 
and ( |1V.2| ), the scattering off a bare nucleon and self-energy like corrections (see, also Fig. 
2, diagrams (a)). These diagrams provide the contribution to the Fermi motion of "dressed" 
nucleons. Besides, there are terms of a pure exchange origin which reflect the fact that 
nucleons in the deuteron are "off- mass-shell" , and terms with renormalization and recoil 
contributions. 



M=±l 



1 



3 P \ (_iR , M n+1 (^)-M n+1 (^ 
2 ° A 2 M n+1 (g?) 



+ 



S3^ 1) (p,p)«S(p)sS 2) (p)*St(p) 



(27T) 



/ 7^Y^m(p)S z VM (-L - iVf (k)) ^(p + k) 



(IV.7) 



(2tt) 6 z n " ' \u(k) 

where S is the total spin of the nucleons, S = \{(T\ + cr 2 ), (with its projection on to the 
z-axis S z ), V n (k) is the one pion exchange potential, (S Z ) D = (1 — 3/2Pd), Pd is the D-wave 
probability in the deuteron, E^(p) = |(E+(p,p) + S^(p,p)) (see eqs. ( |11. 24] ) )- ( [11. 27p ) and 
N®(k), is defined by 



iVf(k) 



1 + 



2m 



1 



N 



2m 



N . 



(IV.8) 



The expressions in eq. ( [IV. 7|) still contain the bare moments, M n (g^). However they may 
be expressed trough the physical moments by making use of 
3 



(l- -P D ) (M n+1 (sf ) - M n+1 (g* 



M„+i(^f; 



= (S z ) D (M n+1 (g? 
(E+ 12) • NU(p)) D (M n+1 (^f) - M n+1 (gi 



N 



(IV.9) 



In obtaining expression ( |1V.9| ), we have used that fact that (M n (g±) — M n (g^)^j ~ g\\ then 



the moments of the deuteron, up to g\ terms, contain only well-defined quantities, viz 



1 / M \ n 

- (-°) x M n+X {g°) = M n+1 (g») AZ 
2\m N / 

+ M n+1 (^)/^3^( P ) (l 



(IV.10) 



Pz 



+ 



rriN 2m 2 
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+ M„ +l(9l »)/g^'(p,k)l 



1 + 



2m jv / \ Ivrifq 



where 



fl A (p) = l E M^(p)(^ + S z + ^- Pz )^(p), (IV.11) 
#*(p,k) = 1 £ A*tf$(p) {<? 2 , K(k)} ^(p + k), (IV.12) 

4 M=±l 

where {S 2 , V n } stands for the anticommutator of S* z and \4 and AZ, which is the difference 
between the renormalization and recoil contributions, has the same origin as in the case of 
nucleons and reads as follows 



P*5?Cp> K(k)| K(v + k). 

(IV.13) 



M=±l 



Figure 2 illustrates the dressing of the moments of the deuteron. The sum of diagrams 
(a) gives the impulse approximation, that is, the scattering off a physical, already dressed 
nucleon. Comparing these diagrams with those in Fig. 1, it is seen that in the impulse 
approximation the helicity amplitude on the deuteron is determined by the scattering am- 
plitudes off free physical nucleons, which ought to be fixed from other experiments. In this 
context, the present approach allows one to avoid the problem as to what amplitudes should 
be associated with "off-mass-shell" nucleon in impulse approximation (see, for instance, Ref. 



42[). In our approach the binding effects are taken into account by terms which are of a pure 
exchange origin (Fig. 2, diagrams (c)) which contain implicitly, via the potential V^(k), the 
contribution of the meson degrees of freedom. The diagrams (b) are the recoil contribution 
terms. 

Applying the inverse Mellin transform to ( [1V.10| ) and omitting the AZ part, the deuteron 
structure function can be obtained in the convolution form 

M D /m 

\£{x)= J f a? [^{tfiv) + (iv.i4) 

where the distribution functions ^^(y) and f^*'(y) are given by 
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&(v) =1 7^3 fh A {v) S (y - 1 - - /r), (iv.is) 



(27r) V m N 2rn N 

d P dk rint, , x 1 



5 l-y+^-)-5 1-y 1 



8(y), (IV. 16) 



where \ IA corresponds to the impulse approximation, or Fermi motion correction, with "on- 
mass-shell" nucleons, and f m * accounts for the binding of the nucleon inside the deuteron. 
Equations ( W.10|) and ( |IV.14|) - ([[V.16|) are the basic result for the determination of the 



moments and the SSF's of the deuteron within the OPE-OBE approach. It will be shown 
later on that eqs. ( |IV.14| )- (|IV.16| ) lead, if proper assumptions are made, to the phenomeno- 



logical convolution model approach used in refs. ||12|| . It is worth recalling here the problem 
as to whether the so-called flux factor has to be considered in the convolution formula for 
polarized deep inelastic scattering |T3| , f43| . In our approach the non-relativistic flux factor 
~ (1 +p z /m N ) comes automatically, as it can be seen from eq. ( |1V.11| ). 

Formulae ( [1V.10D - (|1V.16| ) have been obtained for the pseudoscalar isovector coupling. 



In this case the deuteron wave function ty D (p) appears to be the solution of the Schrodinger 
equation with the one pion exchange iViV-potential Kr(k). Obviously, this wave function 
and the one pion exchange potential are not yet sufficient to describe the properties of the 
deuteron. To this end it is necessary to take into account other mesons contributing to the 
OBE potential, viz. the a, lu, p, rj and 5 mesons [0]. Including these mesons in our approach 
leads to contributions similar to ( |IV.10| ), ( |IV.16| ), except that the wave function \1/ D (p) is 



replaced by the solution of the Schrodinger equation with the full OBE potential. The 
convolution formula can be written in a more compact form by expanding the 5 functions 
in ( |IV.16|) and retaining terms up to kl/m%; one gets 

i AHni ^ KH) .M (IV , 7) 

where 

Md/M 

JA(„\ f dy N " 



y \yJ 



X 



and {S z Vobe)d is the spin- weighted mean value of the potential of the nucleon in the 
polarized deuteron. 
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The second term in eq. (JV7T7) is the correction to the impulse approximation due to 
the binding of nucleons. It can be seen, that this contribution is small (~ ( Vqbe ) l m N) and 
depends on the behavior of the nucleon structure function Qi{x) and its first derivative. 



B. The 3 He 



In this section the Compton helicity amplitude h 



+ M,+ M 



pertaining to the scattering of 



the virtual helicity (+) photon off the polarized (with projection A4 = ±|) 3 He nucleus will 
be presented. Using the Tamm-Dancoff decomposition for the state vector \ 3 He)Q =0t M and 
the operators 0^(n), the following form for the helicity amplitude can be obtained 



h 



2M 



+ M,+ M 



3ff 

m N „ 0.2... 



1 \ n + 1 

E z 



X 



J dV, He ^ e *(p u p 2 ,p 3 ) 



dk gj 1 » 

(27r) 3 4m^a; 2 ^ 



{3kV«(7« (i- - IY«(k)) + (2k z (<r« k) - kV«) t«(7«t« (-L - ^(k))} 



+ (l-Z3 He )E^ 1) (p t )S J + (p,)C'K 



^ e (Pl,P2,P3) 



+ 



(2tt) 3 



' k E f{^-, ^} (- + N$» (k)) + ^?><7« 



dp. 



^ e ( Pi + k, Pfe ,p, -k) 



(IV.19) 
(IV.20) 

(IV.21) 



where 5 2ij = \ (*®C® + a®Cm, dVs He = U 7^(Pi + P2 + P3) and in eq. ([TV72TD i, j 



i=i (27r) f 

list two nucleon pairs exchanging a meson while the remaining fc-th nucleon interacts with the 
incoming lepton, k 7^ In eq. ( |1V.19| ) the reader may easily recognize the corresponding 
"dressing" part for bare nucleons, eq. ( [111.10| ) and ( |1 1 1 . 1 1| ) , which being, included into eq. 
( [LV.20 ), gives the impulse approximation with the physical nucleons. The last term in eq. 



( [IV.21D and the term depending upon the renormalization constant Zs He in eq. ( |IV.20|) 



correspond to the recoil and renormalization contributions, similar to those obtained in the 
deuteron case. The moments of the 3 He SSF's can easily be obtained by recalling eqs. 
( |11.14j ) and (|IH.13| ): the inverse Mellin transform gives the SSF in the convolution form 



He /m N 

r dy 












I y 







X 



(IV.22) 
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where the effective distribution functions f?# e (2/) of nucleons in 3 He contain contribution 
from the Fermi motion of "on-mass-shell" nucleons (impulse approximation) and from the 
nuclear binding 



fi£e(y) = VA(y) + ^.(y)i 



where 



pi 



f/A iv) = / dVs He n?' n (pi, p 2 , ps) 8 ( y - 1 , . . 

ffnl(y) = 



(IV.23) 



(IV.24) 



dV3 He dk 6(y) 

3-<nt.(Pl ; P2, P3, k) 



(27T) 



5(1-3/ + 



2m n 



-5[l-y 



2m 



N 



with the spin-dependent momentum distributions n!j' n and rp^ t being defined by 



n 



f B (Pl. P2. P3) = E A* ^% + (P*) [(1 ± T*(i))/2] * 
A4=±l/2,t 



He 
A4 > 



<At.(Pl>P2,P3,k) 



(IV.25) 



(IV.26) 



(IV.27) 



= E ^ *M 6 *(Pi, P2, Ps) V^S zij [(1 ± r 3 (j))/2] ^ e ( Pi + k, p*, Pi - k). 
Eq. (|TV.22[) can be cast in a more manageable form by expanding the 5 functions as in the 



case of the deuteron 



9[ He {x) 



M 3He /m N 



dy 
y 



gl (lyUy) + 9i [zWiy) 



(IV.28) 



ax m^ dx m^ 

where (a\ ■ Vq BE )3 He is the spin-weighted mean value of the potential of the nucleon ( 

i = p,n) in the polarized 3 He. 

C. Generalization to heavy nuclei 

Let us generalize our approach to heavy nuclei. To this end, let us first discuss the 
unpolarized case. The comparison will be carried up to a given order in (p/mx) 71 , namely 
to n — 2. 
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Guided by the results for A = 2, one can write for a generic isoscalar nucleus A JT8| Q 

dF 2 N (x) (V) A 



F 2 A (x)=Fl A (x)-x 



dx 



m N 



(IV.29) 



where 



-iIA, 



X) 



F 2 N [-\ f A A {y)dy, 



(IV.30) 



and 



fi A (y) 



1 + 



Pz 

m N 



) n(p)6 ly - 



1 + 



+ 



Pz 



2m N m N 



dp. 



(IV.31) 



In Eq. ( |IV.29|) (V)a is the mean potential energy of the nucleon interacting with the lepton; 
it is linked to the nuclear mean potential energy per nucleon (V) = (^a\ J2i<j v ij\^A)/A by 
the relation (V)a = 2(V). The well known relationships between the total energy per nucleon 
€a = Ea/A, the mean nucleon kinetic energy (T) = (^^1 J2iU\^A) /A an d the mean removal 
energy per nucleon (E) are given by €a = (T) + (V) and (E) = 2|e^| + (A — 2)(T)/(A — 1). 

Now we turn to the deuteron eqs. ( |1V.11| ), ( [IVTT^ ) and ( |1V.14| )-( |IVTT6|) . In the definition 
of the distribution functions in impulse approximation, eq. ( |1V.11| ), it is easy to identify 
the non-relativistic analogue of the nucleon spin vector S z + (S • p)/2m^, which may be 
obtained by applying the Lorentz boost operator to the spin vector of a nucleon at rest. In 
the non-relativistic limit the spin vector is the same in any reference frame. Consequently, 
in the extreme non-relativistic limit it should be replaced by S z only. Furthermore, in the 
matrix element of (S ■ p)/rriN in eq. (|I V. 1 1|) only the ^-components of the scalar product 
(S • p) give contribution as it can be checked by direct computation. So that in this case 
the distribution function f IA (p) becomes 

1 



f IA (p) 



M=±l 



m N 



(IV.32) 



3 Note that in eq. (16) of Ref. the "plus" sign in Eq. (16) should be replaced by a "minus 
sign. 
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Expanding the S-f unction in ( |IV.16| ) around the "on-mass- shell" y, S(y — 1 — p 2 /2 



m 



N 



p z /m N ) « 5(1 — y) — 5'(1 — y)(p 2 /2m 2 N + p z /m N ) and substituting in eq. ([[V.16|) the 
difference of two 5-functions by its first derivative we get 



X 



5(1 -y) -8'(1 -y)(p 2 /2m 2 N +p z /m N + E D /m N - (T)/m N ) 



(IV.33) 



In deriving eq. ( [IV.33Q we used the Schrodinger equation to express (S Z ,V) through the 
deuteron binding energy Ed and the kinetic energy (T) of nucleons. It can be easily shown 
that eq. ( |1 V -33| ) can be cast in the form 



dp 



■dEVP(p,E) (1 



Vz 

m N 



Sly 



m N - E - p 2 /2m N + p z 

m N 



(IV.34) 



where we have introduced the deuteron spin dependent spectral function 



7f(p, E) = ^ =1 ( P ) S z ^ =1 ( P ) 5(E - \e D \) 



(IV.35) 



giving the probability to find in the deuteron a nucleon with momentum p, removal energy 
E = \ed\ and spin projection S z . In order to generalize eq. ( |1V.34[ ) to a heavy nucleus, 
we notice that the quantity m N — E — p 2 /2mjv appearing in the 5 function of eq. ( |IV.34|) 
is nothing but the time component of the four momentum of an off-shell nucleon in the 
deuteron, viz. po = Mr, 



p 2 + m 2 N ~ m N - \e d \ - p 2 /2m N . 
Formula (|IV.34|) can therefore be generalized to the case of any nuclear mass number 
A > 2 by substituting the "deuteron spin-dependent spectral function" V^(p, E) with the 
corresponding nuclear spin dependent spectral function |H],fH| . Then eq. ( |IV.34| ) exactly 



coincides with the phenomenological convolution approach used in Ul4",l2"|| (see next Section). 



D. Comparison with the conventional convolution approach 



Using the concept of spin-dependent spectral function [JTJJ], the SSF of 3 He has been 
recently calculated [T2| within the so called convolution approach in which the lepton is 
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assumed to interact with off shell nucleons with four momentum p = (po,p), with po 



Ma — y (E — + Ma) 2 + p 2 , where Ma is the sum of the target nucleus and E = Ma~i + 
mjy — Ma + E A _ X is the nucleon removal energy with Ma-i being the mass of the spectator 
A — 1 system and E A-l its intrinsic excitation energy. Since in our calculations we have used 



the results of Ref. |nj] it is worth comparing the two approaches from a formal point of view. 
To this end, we start with the unpolarized case. Then we can use the analogous of ( |IV.17|) 
for an isoscalar nucleus in the unpolarized case and take its generalization to any value of A, 
as we did in the previous section obtaining eqs. ( |IV.29|) - (|IV.31|) . The convolution formula 



for off shell nucleons reads as follows [13 



F 2 A (x) 



F 2 A - f A (y)dy, 



(IV.36) 



where 



f A (y) 



p 



P(\p\,E)8(y 
m N ) \ 



P 



dpdE 



(IV.37) 



with p + = po + p z and P(\p\,E) being the unpolarized nucleon spectral function. It can 
be seen that whereas in eq. ( [IV.17] ) the binding effect is explicitly displayed through the 
mean potential energy (V), in the conventional convolution approach the binding effects are 
hidden in the definition of the light cone momentum distribution fA{y)- Let us however 
write down the latter in the order p 2 /m 2 N as in Eq. (|IV.i7|) . To this end, we expand the S 



function in Eq. (|IV.37| ) around the point y — 
order p 2 /m 2 N ; we get 



1 + 



+ 



Pz 



2m N TJi jv 



and keep only terms of the 



f A {y) 



1 + 



x 8 [y 



Pz 

m N 
1 + 



n(IPl) 



p- 



+ 



Pz 



+ 



2m N m N 
(E) + (T) A/ (A - 1 



dp 



m N 



S'(y - l) 



(IV.38) 



where n(|p|) = JP(\p\,E)dE and (E) = J EP(\p\, E)dpdE. Placing ( [TVT38D in ([06|) 
and using the relationship (E) + (T) = 2(V) we recover Eq. ( |1V.29| ). Thus we have 
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demonstrated the convolution formulae arising from the OPE and the conventional treatment 
of unpolarized DIS coincide up to the order p 2 /m 2 N \ likewise we have shown that at this 
order the binding effect in the latter approach arises from the average potential energy of 
the nucleon hit by the incoming lepton. Let us turn now to the polarized case and let us 
analyze the 3 He case. The conventional approach yields (TJ| 

f A dy 

J T 

t=n,p 



x 



(IV.39) 



where the light cone momentum distribution G l (y) is given by the following expression: 

Po +P, 



G\y) = JdEjdpPi(p,E)5(y 



(IV.40) 



V m N 

where E = Md + — Ms He + Ed* is the nucleon removal energy {Ed* being the energy 
of the spectator np pair in the continuum), p = Mz He — [(E — m N + Mzue) 2 + |p| 2 ]^ is 
the energy of a bound "off- mass-shell" nucleon, and P\\{p, E) is the spin dependent spectral 
function (cf. eqs. (9) and (16) of Ref. jl2fl). The integral of the spin dependent spectral 
function represents the spin dependent momentum distribution 



n\{P) 



dEPi(p,E) . 



1 + 



+ 



Vz 



(IV.41) 



and considering the 



By expanding the 6 function in ( |1V.41| ) around y 

2mN rriN 

recoil of the two body system non relativistically {Er = p 2 /4mjv), and taking into account 
the differences between the proton and neutron spectral functions (see fT3]|), it can be shown 
that (|TV.22D and (|1V.28|) are recovered. 



V. RESULTS OF CALCULATIONS 



A. The Deuteron 



For explicit numerical calculation of the nuclear SSF, gf (x) and gl He (x), one needs, as 
it can be seen from eqs. ( [1V.14D and ( |1V.28| ), a suitable parameterization of the isoscalar 



nucleon SSF (x). Within the present approach this function contains all the information 
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about the Wilson's coefficients and the influence of the meson cloud on bare nucleon; more- 
over, according to the main assumptions of the effective meson-nucleon theory, it includes 
all the dynamic at distances shorter than the core of OBE-potential. The parameterization 
of g± should be fixed from the experiment; in principle, there exist nowadays experimental 
data of <7i both on the proton [J] and the neutron || structure functions, but they are not 
yet fully complete, especially at very small values of x. In this interval some assumptions 
about the behavior of the nucleon structure function are unavoidable. Moreover, the choice 
of the isoscalar structure function gi(x) determines whether the Bjorken and Ellis- Jaffe 
sum rules will be fulfilled or not. In this sense our results depend on the parameterization 
of the nucleon SSF, however this dependence is found to be insignificant (see below). We 
have chosen here two different parameterizations of the nucleon SSF, describing quite well 
the EMC data on proton and satisfying the Bjorken sum rule. We use the parameterization 



from Ref. [45] as a basic input in all our calculations, and in order to analyse the dependence 



of the results from the chosen parameterization we use the parameterization from Ref. []46 

In our numerical calculations for the deuteron we use the wave function obtained from 
the Bonn interaction, yielding Pjj pa 0.0428 [S3]. Figure 3 displays the numerical estimate of 



the ratio Rd = g± / gf , which illustrates the effects of nuclear structure on g® . The results 
presented in this figure, deserve the following comments: i) within the impulse approximation 



(dotted line; eq. (IV-14) without f mt ), the ratio in the interval 0.2 < x < 0.7 is governed 
by the destructive contributions of the Z)-wave admixture which generates a polarization 
of the deuteron along the ^-direction even if the nucleons have their spins aligned in the 
direction opposite to the polarization. Thus it can be concluded that the effect of the D 
wave is the most relevant nuclear contribution within the impulse approximation; ii) as for 
the interaction term, according to ( |LV.17|) , the binding corrections to the deuteron SSF is 
governed by (S z , Vobeid- Let us first of all estimate the contribution due to pions; a direct 
numerical computation gives 

(S z ,V n ) D ~-5MeV (V.l) 
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(if model ambiguities given e.g. by different choices for the form-factors and wave functions 
are considered, one gets (S z , K-)d ~ —(3 -r- 5) MeV). It can be seen from Fig. 3 that the 
effects from binding due to 7r mesons is similar to the one occurring in the unpolarized 
case, but significantly smaller. It can also be seen that the results are not sensitive to the 
parameterization of the nucleon SSF (cf. long and short dashed waves). We turn now to the 
estimate of the contributions due to other mesons, although it could be anticipated that the 
pions give the dominant contribution to the binding effects, since the pion contribution in the 



deuteron is the most significant one |37|j48|] . In order to estimate the most general boson 
contribution to the deuteron SSF one should, in principle, calculate the matrix elements 
(S z , Vobe)d for all kinds of bosons considered by the Bonn potential model. To this end 
one can use the Schrodinger equation to get 

(S»Vobe) = (l - \Pn) en - (<T) - \{T)^j , (V.2) 

where (T) 0j 2 are the mean values of the nucleon kinetic energy in the S— and D— waves (for 
the deuteron wave function of the Bonn potential we have (T) ~ 10.2 MeV and (T) 2 ~ 4.4 
MeV). Taking into account Fermi-motion and the full interaction effects by eq. ( |V.2| ) results 
in a EMC-like effect, as in the unpolarized case (full curve in Fig. 3). By comparing the 
dotted curve (only Fermi motion) with the full one (Fermi motion plus binding), it can 
be concluded that the main nuclear structure effect in the deuteron SSF comes from the 
presence of the .D-wave, with the binding effect remaining a rather small correction. Note 
that all curves in the Fig. 3 are not shown at values of x smaller then x ~ 0.2. The reason 
is that realistic parameterizations of the nucleon SSF have nodes at small x, which lead 
to "poles" in the ratio / simulating nuclear effects. At values of x smaller then these 
poles, all curves in Fig. 3 tend to the limit (l — |-Pd) as x —> (this is a general result for 
all kinds of parameterization no more singular at the origin than 1/x). 

Figures 4a and 4b display the calculations of the absolute values of the deuteron structure 
function and the comparison with the SMC experimental data: a good agreement between 
our results and the experimental data can be observed. The numerical estimate of the first 
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moment of gf{x) within our approach, J dxg^{x) ~ 0.03, is also in an agreement with the 
experimental result J dxg? {SMC) (x) = 0.023 ± 0.02 ± 0.015 0. 



B. 3 He 



As is well known |L0| the interest in polarized 3 He targets stems from the fact that such 
a system can be considered to a large extent as an effective polarized neutron target. As a 
matter of fact, the two protons in 3 He are mostly in the singlet 1 S'o configuration so that 



the polarization of the 3 He is mainly determined by the polarization of the neutron |TT],|10 
For such a reason DIS experiments off polarized 3 He are aimed at obtaining information on 
the neutron SSF's. However, the non vanishing proton contribution to the total polarization 
and nuclear structure effects could in principle hinder the direct extraction of the neutron 
SSF. We have calculated the 3 He spin structure function g\ He (eq. ( |lV.22j )) using the same 
approximation as in the case of 2 H, i.e. by using eq. ( |1V.28| ). 



The interaction term, obtained using the Schrodinger equation and a three-body wave 
function containing S, S' and D waves, reads as follows 

K • Vg BE ) 3He = 2 {Ip^ - 2 -P D ^ ■ e iHe 
2 2 



2(-(T )sl --(T )D ) 



2pO-E3 He -2{T)\ (V.3) 



e 



2((T) 5 + ^(T) S ,-(T) D ) 

2p n -es He -2(T)? { (V.4) 



where P» = J n||(p)(ip/(27r) 3 , (T)l = J n^(p)p 2 /2rriN dp/(2ir) 3 and €3 He is the mean value 
of the binding energy per nucleon in 3 He. 

The spin-dependent momentum distributions, the effective nuclear polarizations and the 



mean values of the kinetic energies have been taken from Ref. [12], where these quantities 
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have been calculated from the spin- dependent spectral function obtained from the Reid Soft 
Core interaction. The results are 



{< ■ Vg BE )* He « 2AMeV, (V.5) 

W ■ V8 BE )*He « -n-SMeV. (V.6) 
In Fig. 5 the ratio R^He = 9i He I ' 9\ calculated using the proton and neutron spin structure 



functions from Ref. |45] and |^6| is presented, whereas in Fig. 6 the SSF g l He is compared 
with the free neutron SSF. It can be seen that, as in the deuteron case, the contribution 
of binding effects is rather small. It should be pointed out that the results presented in 
Figs. 5 and 6 can hardly be distinguished from the ones obtained within the conventional 
convolution approach of Ref. |L2j (eq. |IV.39|) . According to the conclusions of Section 



4.3 this means that relativistic effects (terms of order larger than ^-) are small, as also 
demonstrated in Ref. [[HJ. Our results fully confirm what found in [|1^], namely that nuclear 
structure effects in DIS of polarized electrons off polarized 3 He are those due to the effective 
proton and neutron polarizations generated by the S' and D waves of 3 He wave functions, 
so that the relation 

gl He (x) « 2p p gl(x)+p n g?(x), (V.7) 

p p = —0.030 and p n = 0.88 being the effective nucleon polarization, represents a reliable 
approximation of eq. ( |IV.39|) at x < 0.6. The smallness of the difference between the 
free neutron structure function g™{x) and g\ He {x) is due to the smallness of the nuclear 
structure effects and is largely independent of the form of the chosen parameterization for 
the nucleon SSF. For instance, using the parameters from the Ref. ||45|| , the results presented 



in Fig. 5 and 6 change by 20%, i.e. by a quantity well below the experimental errors, hence, 
eq. ( |IV.11| ) may be considered as a good approximation for the extraction of the neutron 
SSF. 
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VI. CONCLUDING REMARKS 



The necessity of plausible and precise data on the neutron spin-dependent structure func- 
tion is obvious. Since the information about the internal neutron structure is predominantly 
obtained from nuclear data, usually from the polarized deuteron and 3 He, an appropriate 
nuclear model for subtracting of the effects of nuclear structure is requested. 

In this paper a theoretical approach for the analysis of polarized deep inelastic scattering 
off light nuclei was proposed, which allows a self consistent consideration of the role of the 
Fermi motion and the meson degrees of freedom. Since our model relays on the operator 
product expansion method within OBE approximation, the n— th moment of nuclear SSF 
have been found as a product of two n— dependent functions. The contribution of the impulse 
approximation and of the nuclear corrections to the moments have been separated in an 
explicit form. As a consequence, the inverse Mellin transform yields back the nuclear SSF in a 
convolution form with two distribution functions: one of them describes the electromagnetic 
interaction of the lepton with an on-shell nucleon, whereas the other one describes the strong 
interaction of the hit nucleon with the other nucleons of the nucleus. 

A generalization of the model to the case of heavy nuclei has been proposed and the 
comparison with the convolution approach based upon lepton scattering off bound off-shell 
nucleons |0| was performed with the result that up to order p 2 /m 2 N the two approaches 
coincide. The numerical calculations of the spin-dependent structure functions for polarized 
deuteron and 3 He show that the nuclear corrections are relatively small and essentially 
depend on the spin-orbital structure of the corresponding nucleus, in agreement with the 



finding of Refs. |14| , based upon the conventional convolution approach. For the deuteron 
the main effect of the nuclear structure is due to the destructive role of the orbital motion of 
nucleons with L = 2 and is of the order of magnitude ~ (1 — 3/2Pd). The comparison with 
the SMC experimental results || shows a reasonable agreement of our calculations with the 
data. In case of polarized 3 He the main nuclear structure corrections come from the S'— 
and .D-waves admixtures to the ground state wave function. Since they lead to a partial 
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depolarization of the neutron inside polarized 3 if e and to an effective proton polarization, 
g\ He slightly differs from the free neutron spin-dependent structure function g™. The binding 
effects in both cases are found to be small. These results agree with those obtained within 
the convolution approach [|T4|,|1^.[J . 



To sum up, we can conclude that the convolution approaches so far proposed describe 
fairly well the peculiarities of polarized DIS off polarized nuclei, so that nuclear corrections 
can be estimated in a reliable way. In closing, we would like to point out that our method 
based upon the OPE within the effective meson nucleon theory, allows one to microscopically 
understand the origin of the binding effect which is present in the convolution approach. 
Moreover, since the expressions for moments have been obtained explicitly, this allows one 
to estimate the influence of nuclear effects on the Bjorken sum rule via direct calculations 
of the first moments. 
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FIGURES 

FIG. 1. The dressing diagrams for the helicity amplitude corresponding to the forward Compton 
scattering off the proton. Black dots denote the scattering off a bare nucleon and the dashed 
line the pion propagator. The first diagram is the contribution of the impulse approximation 
and renormalization terms, the second and diagrams represent the recoil effects and the last two 
diagrams are the self-energy-like interaction terms. 

FIG. 2. The helicity amplitude for the forward Compton scattering off the deuteron. The 
notations are the same as in Figure 1. The first three diagrams represent the impulse approxima- 
tion with physical nucleons plus renormalization effects and the remaining two diagrams are the 
deuteron recoil and interaction terms, respectively. 

FIG. 3. The ratio of the deuteron and the isoscalar nucleon spin-dependent structure functions 
(SSF). Curve 1 : the impulse approximation; curve 2: impulse approximation plus interaction term 
{{S z , Ftt})- Both curves have been computed using the parameterization of the nucleon SSF from 
Ref. pip; Curve 3: the same as in curve 2, with the parameterization of the nucleon SSF from Ref. 



[46]; Curve 4: the same as in curve 2, plus the contribution from the interaction term {{S z , Vobe}) 
representing the total contribution af all mesons considered in the One Boson Exchange (OBE) 
interaction. The arrow indicates the value (1 — f-Po) (see text). 

FIG. 4. (a) The weighted deuteron spin-dependent structure function xgi(x); (b) the first 
moment Mi( 5 f ) = f™n/m N g D^ dx 

versus th.6 low6r limit of intogra/tion x m i n . The dotted (full) 
line corresponds to calculations with the parameterization of the nucleon SSF from Ref. |45| (Ref. 
). Experimental data from Ref. ||. 



FIG. 5. The ratio of the 3 H e and free neutron SSF calculated within the impulse approxima- 
tion (dotted) and within the impulse approximation plus the interaction term {{S z , Vobe}) (full). 
Curves 1 (2): parameterization of the nucleon SSF from Ref. (Ref. |j46f). 
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FIG. 6. The SSF for 3 He calculated using the parameterization of the neutron SSF from Ref. 



[45], given by the dotted line. Dashed line: the impulse approximation; full line: the contribution 
of the impulse approximation plus binding effects. 
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